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ABSTRACT 


A theoretical investigation of the effects of short period 
angle of attack oscillation on air-to-ground weapons dive 
delivery accuracy was conducted. The investigation included 
the computer modeling of the longitudinal control system of 
the F-4 aircraft. This longitudinal control system was used 
to fly the aircraft through its standard conventional dive 
bomb delivery maneuver, making corrections for various dive 
angle errors, detected four thousand feet above release 
altitude. 

The dynamics of the F-4 aircraft motion was represented 
by a matrix of dimensional, body axis longitudinal stability 
derivatives and force coefficients. This matrix representa- 
tion of the aircraft dynamics was used to solve for the time 
History Of aircraft. perturbation, velocity, angle of attack, 
pitch angle, pitch rate and flight path angle, following 
correction to the desired delivery maneuver flight path, 

The perturbation angle of attack is a measure of the 
gunsight aimpoint error at any given time, if all other de- 
livery parameters are met, and converts directly to weapon 


impact error. 
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DEFINITIONS OF SYMBOLS 


REFERENCE QUANTITIES 


FRL 


total aircraft velocity in reference state 


aireraftt velocity parallel to FRL in 
reference state 


aircraft velocity perpendicular to FRL 
in reference state 


aircraft reference Mach number 


reference flight path angle 
reference pitch angle 
reference angle of attack 
aircraft mass 

aircraft weight 

air density 

wing area 

wing mean aerodynamics chord 
dynamic pressure 

acceleration of gravity 

moment of inertia about y axis 
thrust incidence with FRL + oh 


perpendicular distance to thrust line 
Lrom, ¢. 2 . 


longitudinal reference axis-fuselage 
xreference line 


longitudinal distance from c.g. to 
bobweight 


(ft/sec) 


(ft/sec) 


(ft/sec) 


(radians) 
(radians) 
(radians) 
(slugs) 
(pounds) 
(slugs/ft?) 
(ft?) 

(ft) 
(1bs/£¢7) 
(ft/sec™) 
(slug- ft?) 


(deg) 


(ft) 


(ft) 





fensithar) 
(Phe iby) 
(#pele) 
teianon) 


(© seh sact«) 


(? 7%) 
(3?) 
Cras 
{*s0e yT} 
(* 9) - pale 
3; 


















Asqits <n tai a 
Ze 8 iit reaaeeh 

sagiow, BE a it rs 

dbs irah Br a 4 

AYR ghiw 

broths 2a is nele be imam yates : 
i cieestige Samson om 


(tivatg. to MOPIETOLGISE = 


xe  tuods 6 besa Po Pasngm: =. 
AHR hw worth ogre pets a 


in yk aes 


stad’ Ia bibdate 
edi) S2neteiag 


PERTURBATION QUANTITIES 


perturbation flight path angle 
perturbation pitch angle 


perturbation angle of attack 


perturbation velocity component along FRL 


perturbation velocity component perpen- 
dicular to FRL 


perturbation pitch rate 


perturbation stabilator position 


TOTAL QUANTITIES 


ALPHA 


M 


tetal slight path angle 
Daa ee ¥ = 2 = ALPHA 
total pitch angle ¢G = 0, +8 


total FRL angle of attack 
ALPHA = a, + a 


total velocity compert along FRL 
= hy tu 


total Mach number 


(rad or deg) 
(rad) 
(rad/milli- 
radians 
when refer 
to sight 
angle) 


(ft/sec) 


(ft/sec) 
(radians/sec) 


(rad-pos, 
Ly, Es up) 


(rad/deg) 
(rad) 


(rad/mils) 


(ft/sec) 


STABILITY AXIS DIMENSIONLESS COEFFICIENTS AND DERIVATIVES 


i) 


D/qS 
dC) /da 


dC, / 3M 


8Cp/356 
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c. = L/qs 

Cho = 9C, / 3a 
Cis = [2V_, /c] 3C,/3a 
C = 3C,/aM 

Ly L 

Gi = [2V_ /e] 3c, /3q 
Lae = 9C,/95o 

Ch = m/qSc 
Ce = aC /3a 

Cs = [2V_ (/¢] 8c / oa 
CaM = aC / aM 

Cua = [2Vp /c]ac/3q 
“ns = 9C/ 35. 


BODY AXIS DIMENSIONLESS COEFFICIENTS AND DERIVATIVES 


Cy = Cy EOSsO, = Cy sina, 

Cy a = (Cy y = Cy) COSHo - (Cp + Cr Sind 
Cyn = Com COSAo - Com sindy 

C == COSdp - C Sind 

Xb, Dé, Lég 

Cy = Cr COSM,) + Cysinado 

Ove = (Cry + Cy) cosag *+ (Ch y = Cy) sind 
& = ¢ Ecosay ot GC sind 

Néo Lé, Dé, 

CNs = Cy 3Sindo 
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Ci = CiMcOS% 9 + Comsineg 


All moment coefficients and derivatives are the same for body 


axis as for stability axis. 


BODY AXIS DIMENSIONAL DERIVATIVES 


























pSUp yg W, 
X ei x 20, Cy) (1/sec) 
pSUG 2W, > 
2 3 M 2 
Xo jot Gi (Cy, Te: + 5 Cyy)) (ft/sec“rad) 
pSU c 
Xa i ia % (ft/sec rad) 
q 
- pSV5. : 
Xs . PAG (ft/sec rad) 
S 2m XS, 
eSU W 
M 0 
ar Abie) es mae 2U Cy ) (1/sec) 
pSUG 2W, 
= - - M 2 
zy ea er eNy “Ue (Cy + 5 Cy) (ft/sec” rad) 
pSU,c 
26 z am = °NG (ft/sec rad) 
pSU oc 
24 io ia Ng (£¢/sec rad) 
eSv2 
2 To 7 
PSq- perpurim:iniN&g (ft/sec? rad) 
pSU 9c M W, 
pSU oc 2W, e ; 
Moy FemEee  ne: Om 2 St? (1/sec“rad) 
pSc 7U? 
MG ~ 41 Vr One (1/sec) 
a 


11 








u 


T 

: C 
41 mq 
pScv., ? ; 

0 

C 
ae mS ¢ 

aa Cd 

mV > oM 


CONTROL SYSTEM DIAGRAM TERMS 


- palot~stick, force 


bobweight stick force 
Fop - FR 
control stick movement 


stabilator movement due to stick travel 


stabilator movement due to stab. aug. sys. 


6 + 6 


SST SSAS 


normal acceleration at bobweight 
perturbation pitch angle commanded 


Oo. =i a(6, 


perturbation dive angle commanded 
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I. INTRODUCTION 


The theory and practice of fixed gunsight, air-to-ground 
weaponry is well established and differs only slightly from 
one aircraft to another due to minor handling quality and 
weapons release system differences among aircraft. Just as 
well established is the suatyeds of impact error sensitivi- 
ties to variation of actual from desired release parameters. 
Table I shows the desired release parameters and error sen- 
Sitivities in the most widely used F-4 dive bomb delivery. 
It is not the purpose of this study to disagree with any of 
the release parameters for this or any other delivery maneu- 
ver, but rather to demonstrate that increased emphasis in 
training on one of the error sensitivity categories, that of 
gunsight placement as influenced by FRL angle of attack per- 
turbations, can greatly increase pilot competence in air-to- 
ground weapons delivery. 

A lack of emphasis on angle of attack input to gunsight 
picture seems justified on the surface, since everything 
which is an input to sight angle can be computed and set 
prior to the delivery maneuver and should not be subject to 
Eicon. (In addition, Table I.shows that the impact error 
due to a sight picture error on one mil (0.1% of a radian) 
is smaller than the impact error due to an airspeed error 
Meio kHOtS, Or a flight path dive angle error .of one de- 


gree. For these reasons and since only one FRL angle of 
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Table I. 


Desired Delivery Parameters and Error 
Sensitivities for Maneuver Studied 


Release Parameters 


Dive Angle 40 Degrees 
Release Airspeed + S00 KTAS 
Release Height 5000 Feet AGL 
Sight Angle 116 Mils 
Release Horiz Range 4720 Feet 
Release Slant Range 6875 Feet 
Parameter Error. Impact Error 

Dive Angle 

1° Steep ‘33 Feet long 

1° Shallow 35 Feet short 
Release Height 

100 Feet High 13, Peet short 

100 Feet Low 1S Feet" tony 
Release Airspeed 

IG KIS Fast So Feet fone 

10 KTS Slow 37 Feet short 
Sight Angle/Aim Point 

= Mii/i Mil Long 10 Feet long 

- 1 Mil/1 Mil Short 10 Feet short 


(thas information was calculated from data in Ref. 3. Sight 
angle includes weapon trajectory drop, FRL angle of attack 


and gunsight parallax correction.) 


attack is associated with one steady flight condition at a 
given weight and configuration, it has been accepted that the 
angle of attack input to the sight setting is a quantity which 
varies only with release airspeed and dive angle. This is 
tinue, flowever, only if the aircraft is free from all longi- 
tudinal oscillations rising from the most Peeeut atmospheric 


disturbance or if longitudinal control action were taken. 
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Most of the emphasis given to sight error due to angle of 
attack perturbations in the past can be summed up in the 
Statement: '"'The aircraft must be in one-G flight at release 
in order for pipper placement to be an accurate indication 
of bomb impact point." Very little emphasis has been placed 
on the fact that a finite time interval must expire after 
the most recent Pare ALE) control input before the actual 
FRL angle of attack matches the steady state angle of attack 
fer the desired flight condition. If bomb release occurs 
during this finite interval, any difference between perturbed 
FRL angle of attack and angle of attack computed for the 
Sight setting will manifest itself as a bomb impact error. 
Finding the magnitude and duration of this angle of attack 
perturbation after making various dive angle error correc- 


tions was the purpose of this investigation. 
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it. - PROBLEM "SET-UP 


The standard F-4 visual dive delivery of the MK 82 
LDGP Bomb with conical fins was chosen for the investiga- 
tion. This maneuver involves a roll-in from 12,000 feet at 
400 KTAS on a heading 90° to desired run-in line. Roll-in 
is accomplished by a 135° bank angle descending turn to a 
40° dive along the run-in line. Roll-in is accomplished at 
military power and completed by approximately 10,000 feet. 
Power is reduced to idle at 430 KTAS (at approximately 9,000 
feet), and the aircraft accelerates at idle power in the 40° 
dive to 500 KTAS at the release altitude of 5,000 feet AGL. 

The largest impact error in range is due to dive angle 
error at release. The causes for this error, shown in 
Figure 1) are< 

1. Ina dive angle other than 40°, with the gunsight 
pipper on the target, the range to the target is different 
than that planned for at 5,000 feet, while the bomb range 
from this point is approximately the same as if in the plan- 
ned dive angle. 

2, Ihe FRL “angle of attack is dependent on dive angle, 
but the input to the gunsight for angle of attack is based 
on a 40° dive angle. 

In this study, pilot detection of a dive angle error at 
Pauvurteere (a typical point, to detect an error, shortly 


after roll-out on final run in line) caused the pilot to 
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make an additional dive angle change to correct for the error 
in dive angle. 

There are two different situations which can be considered 
to be the same dive angle error. A five degree shallow error 
exists for either of the following cases: 

1. The pipper placement is correct for 9,000 feet alti- 
tude but the aircraft dive angle, as noted visually outside 
che “arrerattior on’ the’ airtraft: attitude indicator > Ws" 35°" 

2. Dhe aircraft Wsactially in a’ 40°° divebut the’ pipper 
placement at 9,000 feet altitude is 5° short of where it 
should be on the ground at that point. 

In either: case the aircraft is approximately 1,850 feet 
(horizontal range) farther from the target than it should be 
at that altitude in the maneuver, and the pilot must shallow 
the dive still further to re-position the aircraft on the 
proper flight path over the ground, and then re-establish 
the desired 40° dive angle. Since pilots rapidly acquire 
skill at’ rolling into a’ desired dive angle, the normal: dive 
angle error situation is to be in a 40° dive but find the 
pipper aimed too far short or not short enough of the target 
at the 9,000 feet altitude check point. This was the case 


considered in this investigation. 


A. DIFFERENCES BETWEEN ACTUAL MANEUVER AND COMPUTER SIMULATION 
1. Actual Maneuver 
In an actual run the aircraft is constantly increasing 


airspeed from approximately 430 KTAS at 9,000 feet to 500 KTAS 
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at 5,000 feet with changing lift and drag coefficients. Air 
density increases as altitude débwe dex and yp tral target 
elevation is not at sea level. 
2. Computer Simulation 

In the computer simulation it was assumed that re- 
lease airspeed was attained at 9,000 feet and no further 
acceleration was achieved during the runs except as a result 
of perturbations due to longitudinal maneuvering. Since the 
assumed quasi-steady state velocity is the actual desired 
release velocity, results in the area of release altitude 
were very accurate. Air density was assumed constant at 
standard day density for 7,000 feet M.S.L., the mid-altitude 
for the run. Target elevation was assumed to be standard 
sea level. Other assumptions and conventions used were: 

gear bat earth fixed an ‘space. 

b. Atmosphere fixed to earth (no wind), 

GC. Rigid airfrane, 

d. Body axis system along FRL used throughout, 

e. X-Z plane of symmetry, and 


£. Small perturbation angles and velocities. 


B. REPRESENTATION OF AIRFRAME DYNAMICS 
The following body axis force equation system was ob- 


tained from Ref. 1 in the Laplace domain, 
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where notation and definitions are given in Appendix A. The 
above system is represented in the time domain by the follow- 


ing linear differential equations. 
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Substitution of Equation 2 for a into equation 4 yields 
the following time domain matrix differential equation in 
the perturbation quantities u,a, 6.and q. 
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C. COMPUTER REPRESENTATION OF F-4 LONGITUDINAL CONTROL SYSTEM 
Reference 1 contains diagrammatic representation of the 
individual components of the F-4 longitudinal control system 
as a linear control feedback system which can be combined to 
form the complete linear control system shown in Figure 2. 
Either Block Diagram Algebra or the Signal Flow Graph tech- 
nique may be used to reduce the system in Figure 2 to the 
simplified closed loop system shown in matrix equation form 
in Figure 3. Derivation by Block Diagram Algebra techniques 
of the following Laplace domain representation of the con- 


trol system is presented in. detail in Appendix A. 


(1+G2G3G,HiH2(S) - Gy GsH3aHy (S))5_(S) = G1G2G63G,(S)[y.-y(S)]. 


As shown in Appendix A, precise representations in the 
Laplace domain of all the control system transfer functions 
would lead to a 42nd order differential equation, the numer- 
ical solution to which would be impractical for the purposes 
of this study. For this reason a method was sought to sim- 
plify the more complicated transfer functions which were not 
directiy indicative of airframe dynamics: the aircraft feel 
system, defined G,(S); the stabilator actuator, defined G,(S); 
and the vertical acceleration feedback system to the Bob- 
weight, defined H,(S). These simplifications are shown in 
Appendix A. The 8/6. (S) transfer function, defined G;(S), 
although complicated, was not simplified since it represented 
actual airframe dynamic response to stabilator movement and 


its accuracy was considered critical to the investigation of 
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the problem. With the above simplifications, the following 
equation is the Laplace domain representation for the sta- 
bilator position in terms of pilot stick force x and dive 
angle change required You: (See Appendix A for fee 


tion.) 


feet ee aso ae 55.575" + 95.05" + 1.2018 + 0.0224) 5,(S) = 


- 0.01113K, (S° + 3.802S* + 17.99S? + 16.37S? + 1.2015 


+ 0.0224) [y.-y(S)]. 


In the time domain this becomes the following fifth 


order linear differential equation: 


pet?) +. igang. (7) + 352576.) + 95:06. + 1/2016. + 0.02246, = 
S S S S 


S S 


SiG ammbak ny yj!) *'s 802 ty.-4)"") + 17.99¢y,-7) 


petGeemiye- 7) +°1.202(y.-y) + 0.0224 (y_-¥))- 


The state space technique of Reference 5 for representa- 
tion of nth-order systems of linear differential equations, 
in which the forcing function involves derivative terms, was 
used to cast the control system as a set of five linear first 
order differential equations, in general, for the following 


fifth order linear differential equation, 


8 (3) 


Y 7 acy vu any + aoy ¥ aay vi agyY = 


= g, °>) os g.u(4) + gut) + Bou oe Bau + Bou 


with the following definition of terms 
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differential equations. is 
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In the case of the control system, z is defined to be the 


stabilator position 6 and u is the forcing function © aie 
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The subscripts 5 through 9 were used since the control sys- 
tem equations were to be coupled with the four airframe 


equations of motion to form a ninth-order system for numeri- 


cal solution. 


D. AIRFRAME AND CONTROL SYSTEM COUPLED 

The airframe and control system equation were coupled 
through So» defined by ee oe) Using the” fact’ that 
flight path angle [T is defined by I = O-ALPHA and the flight 
path perturbation angle is defined by y = @ - a, the final 
set of nine coupled linear differential equations appears as 


the teedhack (control system of Figure 3. 
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III. PROBLEM SOLUTION AND RESULTS 


A. CASES CONSIDERED 

Two separate initial conditions were treated. jnin the 
first case the pilot discovered at 9,000 feet that he was 
2.2 degrees shallow of his desired flight path. In the 
second case he discovered that he 1.9 degrees steep of the 
desired track. For each dive angle error two alternative 
techniques were used to re-establish the desired 40° track 
over the ground. The first was to shallow or steepen 
further an amount which would allow re-establishment of the 
desired track immediately prior to the release altitude of 
5,000 feet. The second technique was to shallow or steepen 
more than the first to allow re-establishment on the desired 
flight path early enough for the major protion of the air- 
frame short period oscillations to die out prior to the 
‘release altitude of 5,000 feet. Figures 4 through 7 illus- 
trate the«two error cases and. the flight paths used to cor- 


rect for these errors. 


B. METHOD OF AIRCRAFT CONTROL 

In shallowing or steepening the aircraft dive angle, 
the full system of nine equations was used to solve for the 
time response of u, a, 6, q and y to the commanded perturba- 
tion dive angle Xin@ When the correction dive angle or the 


re-established 40° dive angle was approached, the reduced 
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system consisting of only the Chae airframe equations was 
used to solve for the perturbation quantities. This was 
done for two reasons: 

1. The simplified transfer functions of the longitudi- 
nal control system were invalid in the pitch rate range be- 
tween -0.5° per second and 0.5° per second. Obviously a 
constant dive angle or zero pitch rate falls» within this 
range. 

Z. Sanee tt was, the. short, period. angle, of, attack. per- 
turbation which was of interest, the radian frequency for 
this mode of oscillation was such that airframe damping 
rather than pilot control movement is the factor which 
dampens the oscillation. ’ 

The pilot then).flew,the :aircraft,from, the; initial,dive 
angle to the correction dive angle through the control sys- 
tem, anticipated achieving the new desired dive angle based 
on pitch rate and set«the.stabilator te the, reference posi- 


tion (perturbation 6. = 0) for the new dive angle. 


S 
After the aircraft had flown the required time.on the 
adjusted flight path to approach the correct 40° dive angle, 

the controlling sequence was repeated with the pilot re- 


establishing the 40° dive and the airframe dynamics in the 


new reference state maintaining it. 


C. COMPUTER THEORY FOR AIRCRAFT CONTROL 
Reference 6 contains computer programs for computational 


assistance in the study of linear control theory. The 
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programs used were GTRESP, CALCU, RUNGE, TRESP, YDOT, Y8VSX 
and user supplied subroutine RFIND. Program listings and 
explanation are contained in Appendix D. The purpose of the 
programs is to furnish the graphical time response to the 


following system: 


4 xeCEL raat) ls (bh luct)] 
[u(t)] = K(r(t) - {k}" {x(t)}) 
y(t) = {c}! {x(t)}. 


In the case under consideration, {x} is the vector 
of perturbation quantities, [A] is the nine-by-nine or 
four-by-four plant matrix defining the airframe and control 
system, {b} is the control vector, K = Le is the pilot stick 
force per radian of difference between actual and commanded 
dive angle, r is the commanded perturbation dive angle, {k} 
1s the yeetor,of feedback coefficients, and y is the actual 
perturbation dive angle made up of the output vector {c} 


multiplied by the airframe perturbation vector 
Z (h 
Ce Sr hOs-1 yt adidusart,qi ). 


Following is the detailed breakdown of one of the 
situations considered. At time zero the pilot™at 95000 
feet altitude in a 40° dive determined that he was 885 feet 
horizontal range short of the correct position for that al- 
titude in the delivery maneuver. This was equivalent to 
being 2.2° shallow. He elected to correct for this by shal- 
lowing 5°, flying at a 35° dive angle until approaching the 
desired track, steepening back to a 40° dive angle and flying 


at that 40° angle to release. 
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GTRESP was paved with the following data: 

1. the nine-by-nine plant matrix with airframe data 
referenced to a 40° dive angle and the control system, 

zZ. the control vector, 

3. the output vector, 

4. the fendback coefficient vector, 

5 Ky based on required dive angle change of 5°, 

G. Ve equal to +5° in radian measure, and 

1.) Zero initial perturbation vector referenced to 
40° dive conditions. 

GTRESP ‘then used the other subroutines to perform a 
fourth order Runge-Kutta integration to compute and plot the 
time response of the nine perturbation quantities. 

At 0.975 seconds the perturbation dive angle was 
approaching 5° or IT was approaching 35°. The airframe per- 
turbation quantities for time equal to 0.975 seconds were 
transformed to perturbation quantities for a 35° dive ref- 
erence state, and GTRESP was called with the following in- 
put data: 

1. the reduced plant matrix of airframe data only 
referenced to a 35° dive, 

Zs the output; vector, and 

ee abreratt perturbation quantities referenced to 
35° dive at time equal to 0.975 as the initial condition 
vector. 

Failure to input the control vector, feedback coeffi- 


e1ents, K. and Ye bypasses the control system which is not 
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accurate for zero and very small pitch rates and simulates 
the pilot controlling the aircraft ip Seuedang and maintain- 
ing the zero perturbation stabilator position for the 35° 
dive reference condition. Again GTRESP computes and plots 
the time response of the aircraft perturbation quantities. 

At 7.2 seconds total elapsed time the aircraft was 
approaching the desired 40° track over the ground and GTRESP 
was called with the large plant matrix for control back to 
a 40° dive. Input was the same as for the first maneuver 
except that the aircraft dynamics portion of the nine-by- 
Nine matrix was based on the 35° dive since that was the 
condition from which perturbations occurred, and the ini- 
tial condition perturbations remained referenced to the 35° 
dive. 

GTRESP computed and plotted the perturbation quan- 
ebeies. curougn the transition back to a’40° dives At 8.175 
seconds total elapsed time the 40° dive on the correct 
flight path was re-established and it was simulated that the 
pilot neutralized longitudinal control by calling GTRESP, 
with the reduced plant matrix and initial perturbations 
transformed to those for a 40° dive reference state, to 
calculate and plot the perturbation dive angle, velocity, 
angle of attack, pitch angle, and pitch rate to release 
point. ! 

On this particular run it was arranged to have the 
aircraft re-achieve the desired 40° dive angle at exactly 


the 5,000 feet release altitude. Therefore, the angle of 
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attack perturbation, converted to mils, existing at 8.175 
seconds is in fact the release sight picture error that the 
pilot would have had on this maneuver at release altitude. 
Pertinent aircraft position and gunsight error information 
concerning each of the four bomb delivery maneuvers simu- 
Hated 3s Jasted an Table IT. 

Gompeberentput showing time history in tabular form 
of aircraft perturbation during the four phases of each of 
the four runs tested and in graphical form for the final 


delivery phase is enclosed in the Computer Output Section. 


D. RESULTS OF THE FOUR RUNS 

Table II shows that the total time available for ma- 
neuvering after detecting a dive angle error, 4,000 feet 
above release altitude, is very small for a 40° 500-knot 
dive bombing run. Therefore, any corrections made must be 
both timely and correct, if they are to succeed in increas- 
ing the accuracy of the bomb impact. Table II also shows 
that by merely increasing the correction dive angle from 
5° to 8°, and thereby re-establishing the desired dive path 
1.75 seconds prior to release rather than at release for 
the shallow run and 1.4 seconds prior to release rather than 
at release in the steep run, reduces the gunsight position 
error from 83 to 9 mils for the shallow case and from 80 to 
Smits for the steep case. 

Examination of the graphical time histories of perturba- 


tion angle of attack, presented in Appendix C, shows that 
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any additional time spent on the desired dive path rapidly 
decreases the angle of attack perturbation, thereby rapidly 
increasing the accuracy of the gunsight position relative 


to actual bomb impact point. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


Table II shows the elapsed time in each of the four 
weapons delivery maneuvers tested at which the aircraft 
finally becomes established in a 40° dive on the proper 
track over the ground for an accurate bomb hit. The Conm- 
puter Output shows the FRL perturbation angle of attack at 
numerous time increments in each run. It is obvious for 
both the initially steep and initially shallow runs that 
the 5° correction dive angle, although sufficient to correct 
mMicnaitcratt back to desired dive path by release’ altitude, 
Zeriot Large enough to correct’ back to the desired dive 
path soon enough to allow for sight picture settling time, 
due to angle of attack oscillations of the short period 
longitudinal mode. However, by merely increasing slightly 
the magnitude of the correction dive angle change, it was 
‘possible to reposition the aircraft on the desired flight 
path early enough before the release point to allow time for 
the sight picture to stabilize sufficiently to achieve an 
acceptable bomb impact. Examination of the graphs showing 
time history of angle of attack perturbations during the 
final phase of the delivery maneuver shows that even earlier 
establishment on desired dive path is desirable and would 
increase the accuracy of the sight picture further. 

In actual practice, a pilot on detecting a dive angle 


error would not consciously attempt to hold a dive angle 
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correction for a specific amount of time with the idea that 
in so doing he would re-establish his aircraft on the proper 
flight path. If he could do this, there would be no need 
for the gunsight. All the cases described herein would 
actually achieve very accurate bomb impacts, since for an 
aircraft at the correct velocity and in the desired dive 
angle along the proper path over the ground, very large 
angle of attack perturbations produce bomb impact errors 
Only insotar as they change the aircraft dive angle. The 
error which would actually be produced in the case where 
all delivery parameters are met except for an angle of at- 
Lgeck Crror would be that the gunsight reticle would be far 
err the target. 

However, since the pilot in a visual dive bomb delivery 
uses visually acquired data for dive angle and aircraft 
positioning, he is entirely dependent on his gunsight for 
accurate aircraft positioning information; and during a 
period when he is receiving inaccurate gunsight information, 
he has no good indication of his position relative to the 
desired release position. 

It has long been standard practice in visual air-to-air 
gunnery with a gyro-stabilized gunsight, which develops its 
own lead angle based on the normal acceleration of the fir- 
ing aircraft required to track the target, to teach students 
not to squeeze the trigger until they have maintained a 


smooth tracking solution for two seconds on the target. This 
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ee to allow angle of attack oscillations from a changing 
pitch rate to die out prior to firing. 

The need to apply this same technique to visual air-to- 
ground delivery is just as important, or even more so when 
it is considered that a gyro-stabilized sight has inertia 
to decrease the magnitude of the effect’ of aircraft angle 
of attack perturbations; whereas the sight line of a mech- 
anically fixed sight moves directly with angle of attack 
perturbations an amount equal to the movement of the FRL. 
Nor does a pilot in air-to-ground delivery have the option 
to wait until perturbations die out if Sio8 will take him 
below the planned release altitude, since the computed sight 
angle is good for only the planned delivery parameters. 

What he must do is to detect dive angle errors sufficiently 
early in the run and make corrections sufficiently large 
that all corrections are completed and all significant 
perturbations eliminated prior to the planned release height. 

This; study, pertained specifically. to, the, F-4,aircerafe, in 
a 40° dive, MK 82 LDGP Bomb, delivery. However, since all 
Navy fighter and light attack jet aircraft have similar 
short period characteristics with slight variance in damping 
ratio and oscillation frequency, and since all fixed-sight 
air-to-ground weapons delivery maneuvers have similar error 
sensitivities (differing only in the magnitude of the range 
Srror tor each mil of aimpoint error), the results are ap- 
plicable to all fixed gunsight, visual, air-to-ground weapons 


deliveries in all Navy tactical aircraft. As such, it is 
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recommended that similar studies be made to encompass all 
air-to-ground weapons deliveries from Navy tactical aircraft 
and the results be used in the lecture portions of air-to- 
ground training in both fleet and replacement squadrons. In 
addition a graphical display of the results herein can be 
easily demonstrated to a pilot in flight in the following 
manner. Trim the aircraft for steady ‘Tevel flight with the 
gunsight set on the distant horizon. Make a small but sharp 
longitudinal displacement of the control stick and return it 
to neutral. Notice that, although the change in aircraft 
track (speed and altitude) is extremely small, the magnitude 
of the gunsight oscillation above and below the horizon is 
Sipnificant:—and ‘that)-while. the time to damp is short, it 1s 
long enough to lose several hundred feet, if it were in a 
40° 500-knot dive. 

Awareness of this phenomenon can do much to improve the 
proficiency of the not atypical pilot who, on being ques- 
tioned about pipper placement after an unsuccessful air-to- 
ground weapons training flight, is not really able to tell 
an instructor where his pipper was positioned at release on 
most of his runs, this because of late corrections resulting 


in receipt of meaningless data from the gunsight. 
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APPENDIX A: DERIVATION OF LINEAR CONTROL FEEDBACK COMPUTER 
REPRESENTATION OF F-4 LONGITUDINAL CONTROL 
SYSTEM 


The longitudinal control system of the F-4 aircraft is 


shown below in diagrammatical form, 


= 





where, 
4 1 inches stick travel 
Feel System = G,(S)=27.1[ ————————_- GS orca 
2 S44 5.637S+565. 5 PoUiay ores 


Control Stick Gearing = G,(S)=-0.05864 ( radians stabilator travel) 


inch stick travel 


: i z 1 
Stabilator Actuator = G,(S)= 20 ob S420 ] 


Bobweight = H,(S)= 0.1663 ( ounds pepues ht force ) 


GACLAanSaStabakatowe tieavew: 


2 L S 
ee oy Ot I ee) { padian/sec jpitchgate 


With the addition of airframe dynamics and a pilot, modeled 
by KS pounds of stick force per radian of error becween actual 


and desired pitch angle, the representation becomes 
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where the additional transfer functions are, 


Pilot force/radian eg = G,(S) = K (pounds force) 
P radian 


ar/8< transfer function 


= H, (S) = Hees i +1,213S~+27.98S°+2.272S+0.04339 g 


So eideo is aoe tseen azz4 9 oe 
8/6, transfer function 
2 . 
= G.(5). = -24. 64 [—g—2> “4-28 *9 507 ___ pad eat) 
S47) 8006-15. 196-41 1179540, 0224 
/® transfer function = H,(S) = S cradian/sec, 
a Poa radian 


- A block diagram algebra reduction of the control system 


as described in Reference 5 follows: 
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= E,G,G,G, = GoG,G,E, = G,6,G,HjH, 6. 
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Rh tera ~ tment hay 


atest, Raemaiene aes algerie endehonchitimtnbiin 


:  Wikes eae sendeirenen a 
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rea ome 
ae 2°34 


S rod ad G,G,G,HjH, = G,G HH, 


G,G.e7G,(0_ = 6) 


ee £723 
Say eee eis (GG Han, 
(1 + G,G,G,H|H, - G,G,H,H,) 6,(S) = G,G,G,G,[ (0-8) (S)] 


In this final form the control system in the Laplace 
domain can be readily transformed to its state space repre- 
sentation of a linear differential equation for: 6, as a 
function of time and GR =/@).- However reducing! all tie 
transfer functions which are polynomial fractions in S to a 
least common denominator would result in a forty second 
order linear differential equation:~"Theretoere>; representa- 
tion of some of the transfer functions that would yield a 
resulting differential equation more amenable to numerical 


solution but still accurate, was sought. 


a BS8 Ss G, (S) was not disturbed since it represented 
Bet airframe dynamics which was under investigation. 
Neither was the Stab. Aug. System = H,(S). However, a fre- 
quency response investigation of the control stick = G,(S), 
the stabilator actuator = G,(S) and ay = H,(S), using the 
FRESP computer program of Reference 6, yielded Bode Diagrams 


of transfer function magnitudes and phase angles showing 


that over the pitch rate range of interest (0.01 to 0.2) 
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radians/second - a normal acceleration range lupate 2hdt2) the 
above transfer functions could accurately be represented by 
constant magnitudes with zero phase angle. 

With these simplifications, the transfer functions are 


as follows: 


G,(S) = K, H,(S) = 1600 

G,(S) = 0.04792 H,(S) = 0.1663 

G,(S) = 0.05864 H,(S) = S , 

a = 1 
G,(S) = 1.0 H,(S) = 0.15[ sy 1 
g? 4 495 4 5/367 

G.(S) = -24.64 [| ———___5, * 2-25 + 5.567 ______j 

s* + 2.8028” + 15.19S“ + 1.179S + 0.0224 


The Laplace domain representation for b¢ as a function 


of time and (6. - 6) becomes 


4 3 


(See 4ss” + 35.575 2 


+ 95.0S 
4 


+ t. 2008 +0 .0224) 5, (S) = 


3 2 


= 0.01113K, (S° PseoUo mate soon Ti lOnaro © 1.2008 


+ 0.0224)[6. = 01S als 


The value of Ky» pilot force per radian of pitch angle 
.error, was varied depending on what gain was required to 
achieve an average pitch rate of 0.115 radians per second 
which is equivalent to a change in normal acceleration from 
steady state conditions of 3g (typical in fact of the pitch 


rate a pilot would use in making a dive angle correction). 
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In addition, the pilot was assumed to be aware that 
actual flight path angle (dive angle) lagged the pitch at- 
titude angle during a changing pitch rate by the amount of 
instantaneous angle of attack. Since the purpose was dive 
angle change not just pitch attitude change, @6-a or y was 
fed back in the control system to solve for 6. as a function 
of time and LY one? rather than (6-8). THIS “Is the Tepre- 
sentation shown in the matrix feedback differential equation 
of Figure 3. 

The simplification of the transfer functions listed 
above necessitated using airframe dynamics alone, in a 
Simulated stick fixed condition, rather than coupled with 
the control system after the desired dive angle was approach- 
ed, since very small pitch rates were outside the limits for 


which the simplified forms of G(S), G,(S) and H,(S) applied. 
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APPENDIX B: NUMERICAL VALUES OF CONTROL VECTOR AND ELEMENTS 
OF THE PLANT MATRIX FOR REFERENCE DIVE ANGLES 
USED 

The elements of the control vector derived in Sections 
II-B and II-C are invariant with dive angle. Values of 
elements of the control vector are listed in Figure 8. 

Derivation of formulae for the elements of the plant 
MNaerix LS in Sections Ti-=-B\ and il=C, and the matrex 41s 
shown in symbol form in Figure 3. Definitions of the ele- 
Hents. is in “Definitions of Symbols" at the front of this 
study. Calculation of the airframe dynamics entries was 
accomplished in the following manner: 

1. Dimensionless stability axis data were obtained 
from References 1 and 2. 

2. Dimensionless body axis figures were computed from 
the "Definitions of Symbols", using stability axis data 
and reference state angle of attack. 

3. Dimensional body axis figures were calculated using 
dimensionless body axis figures, dimensional derivative 
definitions, and applicable reference state data. Tables 
III and IV are compilations of the unperturbed 40° reference 
state conditions, and aerodynamic coefficients/stability 
derivatives for each reference dive angle used. Figure 8 
lists the values of the plant matrix elements for each dive 


angle. 


49 














aE etotew Odd Bee pepe wiih i a 
‘ole odd to enotriai tee” epic nd wok Lodangs 
Bite 10 ywowt ade fe “Elem te aioe ser a 
2hw segrice eotwanth opaaliate a te aor | 
Tene ghingl (on afd 
 beakesde ove; axeh aa vr Hidesh ‘adet noi ensiiel 

| ®t’ tenon 
newt Losveqmes ovsy aeraged 2tke VbSe Bieta lane 
bieh etxa \Jilidate guise eiedae go ances 7 
fonite F_ oLigiee gtat< oon a 
gnieuv Ooisluoles S234 eettig hehe wee taneteneaie 7 


ovitavitoh Teewiegeelh \ootage? geen Vie eaelanee 


a 


" oe 

ioidal Sib) e338 Sete OY OR eee Bae ener i 
noTgicy Ob bo losovagay sytt to Anite igen oan Wi Bee 
Vi iltdare\ ete Iso. olvatghotam bie. Saotethnes oleate 


(it *oayatoe aoa eet tov 70Ne 


ve.’ 
‘tom Jasly ote 2 eee eae ‘sas 





Since the aircraft was assumed in a constant speed 
(except for perturbations) 500 knot dive, rather than ac- 
celerating from 430 KTAS to 500 KTAS, lift and drag coef- 
These 


ficients from References 1 and 2 were invalid. 


coefficients were calculated for each dive angle from free 


body force balance diagrams. 
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Table III. 40° Dive Reference Conditions and Stability Axis 
Aerodynamics Coefficients and Stability Derivatives 


To . = -40° = -0.6981 radians 
Vr = 844.4 ft/sec 
0 
U = 844.4 ft/sec 
Wy = =0 (58309 ft/sec 
M, = 775 
g., = -39.77° = -0.6941 radians 
a, = 3.99 mils = 0.2286° = 0.00399 radians 
m = 1210 slugs 
W =. 38920 1b 
0 = 0.001933 slugs/f£t> 
S = §30) £7 
=. = 46l04 ft 
q@ = 689.1 1b/tt? 
g = 32.17 ft/sec? 
I, = 122200 slug £7 
BY id eA 
ieees= = 0.57 ft 
ie = a0r59.4 ft 
Ey ¢= 0,06835 
Bie = 201464 
Beato 
Cys. = 0 
Ge = 0.1391 
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Table IV. Body Axis Force Coefficients and Stability Deriva- 


tives. 
t =<S2° 35 -40° -45° -48° 
Cy 0.05026 0.0611 0.06479 - 0.07144 O07 717 
Cy Y -0.002913 -0.002064 -0.006367 -0.006422 -0.01468 
Cyyy 0 0 0 0 0 
Ps. -0.001416 -0.001416 -0.001416 =0.. 001416 -=0), 001426 
Cy 0.1259 0.1517 0.1394 0.1495 0.1628 
Cvy 5.802 SA 22 5.425 3.432 52437 
Cus Oar 0.7 O27 0.7 O07 
Cm 0 0 0 0 0 
CNq Behe Diaeie AE se) TARA Zatz 
CNé. 0.355 | 02555 0.555 0,55 07555 
xX -0.03541 -0.04345 -0.04633 -0.05129 -0.05566 
xy 0.6541 0.4344 ETE) ier ja a 4.534 
xo -0.0254 -0.0254 -0.0254 =0'.0254 -0.0254 
Xb, 0.4274 0.4274 0.4274 0.4274 0.4274 
zy =0.08029 -0.08759 -0.09478 -0.1061 =O), 108 
Ze Pee = 1022. = LOSA, = 1045 = 056% 
Zs =2.007 Ze 7 =Z..007 =2.007 =2.007 
Zo =6.365 = 50/5 =6.2565 =6.3565 -6.365 
Z8q E072 spi Oy ine =DOy v2 =207'2 =107.2 
My -0.001317 -0.001317 -0.001317 =“ 00LSL7. 05001517 
Mi =5.576 -13.76 -13.76 -13.76 =15.70 
M: =0..4553 -0.4553 -0.4553 =-0.4553 -0.4553 
M, -1.079 -1.079 -1.079 -1.079 Se 
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PLANT MATRICES FOR EACH REFERENCE DIVE ANGLE AND THE CONTROL VECTOR 
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APPLICABLE TO ALL DIVE ANGLES 


0.1655 


ELEMENTS OF THE CONTROL VECTOR, 
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FIGURE 8 (CONT.) 
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COMPUTER OUTPUT 


5 DEGREE CORRECTION 
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TIME RESPONSE OF PERTURBATION QUANTITIES DURING PULL UP FROM 40 TO 35 DEGREES 
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OF PERTURBATION QUANTITIES DURING PULL UP FROM 45 TO 40 DEGREES 
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TIME RESPONSE OF PERTURBATION QUANTITIES CURING 40 DEGREE DIVE 
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RESPONSE OF PERTURBATION QUANTITIES DURING PUSH OVER FROM 40 TO 48 DEGREES 
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TIME RESPONSE OF PERTURBATION QUANTITIES CURING 40 DEGREE DIVE 
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